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ABSTRACT 

Type I X-ray bursts are produced by thermonuclear runaways that develop on accreting neu¬ 
tron stars. Once one location ignites, the flame propagates across the surface of the star. Flame 
propagation is fundamental in order to understand burst properties like rise time and burst os¬ 
cillations. Previous work quantified the effects of rotation on the front, showing that the flame 
propagates as a deflagration and that the front strongly resembles a hurricane. However the ef¬ 
fect of magnetic fields was not investigated, despite the fact that magnetic fields strong enough 
to have an effect on the propagating flame are expected to be present on many bursters. In this 
paper we show how the coupling between fluid layers introduced by an initially vertical mag¬ 
netic field plays a decisive role in determining the character of the fronts that are responsible 
for the Type I bursts. In particular, on a star spinning at 450 Hz (typical among the bursters) 
we test seed magnetic fields of 10^ — 10^° G and find that for the medium fields the mag¬ 
netic stresses that develop during the burst can speed up the velocity of the burning front, 
bringing the simulated burst rise time close to the observed values. By contrast, in a magnetic 
slow rotator like IGR J17480-2446, spinning at 11 Hz, a seed field > 10® G is required to 
allow localized ignition and the magnetic field plays an integral role in generating the burst 
oscillations observed during the bursts. 

Key words: MHD - methods: analytical - methods: numerical - stars: neutron - X-rays: 
bursts - X-rays: individual: IGR J17480-2446. 


1 INTRODUCTION 

The Type I Bursts are X-ray flashes that last for tens to hundreds 
of seconds and are extremely bright, being able to produce up to 
10^® erg/s. They are observed in approximately 100 sources, which 
are all low mass X-ray binaries where the accreting compact object 
is a neutron star (NS). The bursts are thermonuclear runaways that 
ignite in the accreted layers on the surface of the neutron stars (see 
Strohmayer & Bildsten 2006; Galloway et al. 2008, for reviews). 
It is usually thought that ignition takes place at one position and 
then propagates across the surface (Shara 1982). Thus, the observed 
lightcurve is the result of the convolution of the local emission evo¬ 
lution at every point on the surface and the propagation of the flame. 

Type I burst provide a very promising means to extract infor¬ 
mation about the underlying NS, such as mass and radius that could 
put constraints on the equation of state of the matter in the core 
of the star (Strohmayer & Bildsten 2006; Miller 2013). However, 
many questions still remain open about the bursts, for example re¬ 
garding what the most likely ignition site is, what nuclear reactions 
take place and how does the flame spreads over the surface. In or¬ 
der to disentangle the information we seek from the lightcurve, it 
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is fundamental to understand how the flame propagates. The prop¬ 
agation mechanism is the focus of this paper. 

Early discussion about possible mechanisms and the result¬ 
ing flame velocity can be found in Wallace & Woosley (1981), 
Fryxell & Woosley (1982a), Fryxell & Woosley (1982b), Bild¬ 
sten (1995) and Zingale et al. (2001) (see also Simonenko et al. 
2012a,b). Spitkovsky et al. (2002) (from now on SLU02) suggested 
that localized ignition relies on confinement of the hot fluid. If not 
confined, the hot fluid would spread out over the surface, dissipate 
its heat and as a consequence no ignition would take place. In the 
presence of confinement, on the other hand, the heat would not dis¬ 
sipate and the flame would spread thanks to some heat transport 
mechanism at work between the hot and the cold fluid. Since most 
of the bursters spin faster than ~ 100 Hz, the authors suggested that 
the Coriolis force could provide the necessary confinement creat¬ 
ing hurricane-like systems. The horizontal scale of both the ignition 
and the propagating front is that of a Rossby deformation radius, 
Rr = s/gH/2Q,, where g is the gravitational acceleration, H the 
scale height of the fluid and H the angular velocity of the star. This 
was convincingly demonstrated by ab initio numerical simulations 
of Cavecchi et al. (2013, 2015). SLU02 also discussed the frictional 
coupling between the top and bottom layers and this is very impor¬ 
tant for this work (see below), but limited themselves to parametriz- 
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ing its effect on the flame velocity. The simulations of Cavecchi 
et al. (2013, 2015) obtained propagation time scales in good agree¬ 
ment but generally somewhat longer than the observed rise times. 
Furthermore, they found that the principal mechanism responsible 
for the heat transport was the top-to-bottom conduction across the 
flame front. These simulations were performed without any mag¬ 
netic field, which we show in this paper dramatically changes the 
velocity and the character of the flame. 

Moreover, Cavecchi et al. (2011) discussed a very peculiar 
source, IGR J17480-2446 (from now on IGR J17480, see also 
Linares et al., 2011, Motta et al., 2011), which is spinning at ap¬ 
proximately 11 Hz. This source is very important, since at 11 Hz the 
Coriolis force is not effective in confining the hot fluid, the Rossby 
radius being greater than the star’s radius. Nonetheless, IGR J17480 
shows bursts and burst oscillations. The latter are fluctuations in the 
burst luminosity that are caused by an asymmetry in the emitting 
pattern modulated by the rotation of the star (first discovered by 
Strohmayer et al. 1996, they are seen in many bursters at different 
rotation rates, see Watts 2012 for a review). IGR J17480 is an ac¬ 
cretion powered millisecond pulsar: fluctuations are also seen in the 
persistent emission, which implies that the accretion flow is chan¬ 
nelled by a significant magnetic field (see Patruno & Watts 2012, 
for a review on the behaviour of accretion powered X-ray pulsars). 
Finally, the burst oscillations and accretion powered oscillations of 
IGR J17480 are in phase. Based on this evidence, Cavecchi et al. 
(2011) proposed that the magnetic field could be the source of flame 
confinement in this case and estimated that a field > 10® G was 
necessary. 

In this paper we explore the effects of the magnetic field on 
the flame propagation of Type I Bursts in order to establish its in¬ 
teraction with the Coriolis force. We explore conditions that apply 
to the bursting X-ray pulsars, where the extra coupling between the 
layers provided by the magnetic field threading them is likely to 
affect the speed of flame propagation, increasing it and thus bring¬ 
ing our results to better agreement with the observations. We also 
dedicate a section specifically to the case of slow rotators like IGR 
J17480 where the magnetic field is needed to provide the extra con¬ 
finement that the Coriolis force is not able to provide. In Section 2 
we present the results of analytical calculation and the numerical 
simulations and in Section 3 we draw our conclusions. 


2 THE INTERACTION BETWEEN THE MAGNETIC 
EIELD AND THE BURNING FRONT 

Here we will describe the results of our analytical calculations and 
numerical simulations, but first we summarize the numerical setup 
common to all runs. 

2.1 Numerical setup 

We use the code described in Braithwaite & Cavecchi (2012) and 
Cavecchi et al. (2013). An important assumption behind our sim¬ 
ulations is hydrostatic equilibrium. This approach is particularly 
well suited for fluids whose horizontal stretch is much greater than 
the vertical extent, as in the case of the accreted layers of NSs. 
We use the so called a coordinate system (see Kasahara 1974), a 
kind of pressure coordinate system. In the a coordinate system, 
the physical height of the fluid is fixed at the bottom to a con¬ 
stant value (0 cm in our case), while for every other grid point the 
physical height will change according to the fluid evolution. On the 
other hand, the pressure at the top is held constant throughout the 


simulation, while the bottom pressure will change according to the 
motion of the fluid. Of the various coordinate systems used in at¬ 
mospheric physics, this is the logical choice in this context, as the 
computational domain can expand and contract to arbitrary degree 
in the vertical direction in response to the heating and cooling of 
the ocean.® 

The code uses finite differences and an explicit third-order 
low-storage Runge-Kutta scheme (Williamson 1980) as the time 
stepper. The time step is limited by a Courant condition (see Braith¬ 
waite & Cavecchi 2012). In our simulations the most limiting fac¬ 
tors to the timestep are fluid motion and the magnetic field evolu¬ 
tion. Spatial derivatives and integrals are fifth order accurate, while 
interpolations are sixth order (for details see Lele 1992 and Braith¬ 
waite & Cavecchi 2012). Variables are evaluated on a staggered 
grid: thermodynamic variables are defined at the centers of each 
grid cell, while velocities and magnetic field components are face 
centred. The staggered grid ensures good conservation of mass, 
energy, momentum and divergence of magnetic field (Braithwaite 
& Cavecchi 2012, for the simulations in this paper, the average 
VB/Bo is at maximum 10“®®). The relation between density, 
pressure, temperature and composition is calculated with the rou¬ 
tine helmeos® of Timmes & Swesty (2000) where electrons can 
be arbitrarily degenerate and relativistic, ions are considered ideal 
gas and radiation pressure is also taken into account. 

In this paper we use Cartesian coordinates, as in Cavecchi 
et al. (2013), but we set our coordinates system such that x is point¬ 
ing in the vertical direction, z in the horizontal direction parallel to 
plane of the simulations (the direction of flame propagation) and y 
is pointing out of the plane. We also neglect one horizontal direc¬ 
tion, assuming symmetry along the y direction. This configuration 
will help us elucidate most clearly the role of magnetic field in set¬ 
ting the fronts propagation speed. Vertically, our fluid is initially 
between Pt = 10®® erg cm“® and Pb = e® *® x 10®® erg cm“®, 
where Pr and Pb are the pressure at the top and at the bottom of the 
simulation. The value of 1.89 in the exponent for Pb implies that 
we are simulating approximately 1.89 scale heights. Horizontally, 
the domain is 6 x 10® cm wide. 

As in previous works, we set up an initial temperature profile 
which is meant to trigger the burst. It is independent of the vertical 
coordinate, but it has a horizontal dependence which is: 

T = To + , ^ 1 (1) 

1 + exp[(z — zo)/dz\ 

The parameters we use in this paper are To = 2 x 10^ K, 5T = 
2.81 X 10^ K, zq = 9 X cm and Sz = 3.6 x 10^ cm. zq 
corresponds to the position where the temperature is To + 6Tj2, 
while 6z is approximately the distance over which the temperature 
transitions from its maximum To + 6T to its minimum To. 

We only consider helium burning via the triple-a reaction 


^ Due to the change of coordinates, the equations for the magnetic part 
described in Braithwaite & Cavecchi (2012) should contain terms propor¬ 
tional to powers of the horizontal derivatives of the height (the terms dh(l> 
in Braithwaite & Cavecchi 2012). We neglect these since they turn out to be 
small due to the extreme aspect ratio of our simulations. 

^ Publicly available at http://cococubed.asu.edu/code_pages/eos.shtml. 
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(Gumming & Bildsten 2000)® 

Q„ = 5.3 X 10®®p5 gjg g-i ^ 2 ) 

and we simulate a fluid which is a mixture of helium and carbon. 

We assume that the bottom boundary is not moving (no-slip 
condition). This is based on the consideration that the density of the 
burning layer becomes much smaller than the density of the ashes 
on which the layer rests. Therefore we allow for a buffer zone of 
pure carbon at the bottom of the simulation which is not affected by 
the flame propagation. In order to do so the initial composition is 
constant horizontally, but the helium fraction decreases downwards 
according to 


1 -b exp[-(a - 0.8)/1.6 x IQ-^] 

(ct = 0 at the top and cr = 1 at the bottom). Erring on the side 
of caution, we also impose a damping force to prevent motion in 
the carbon layer near the boundary. Reassuringly, in later tests it 
turned out that switching off the damping force makes very little 
difference'^. We calculate the acceleration terms dtU^. and dtUy 
and then multiply them by 




1 cr < 0.82 

cos^ f) 0.82 < cr < 0.92 

o' cr > 0.92 


(4) 


The choice of this functional form for x is made to ensure smooth¬ 
ness and the fact that the damping only acts within the buffer zone 
of carbon. 

In our initial configuration the magnetic field is purely vertical 
and has a homogeneous value Bo, which is different for different 
simulations. The gravitational acceleration is p = 2 x 10®'* cm 
s“® and we use the plane-parallel approximation and a constant 
Coriolis parameter / = 212 (the /-plane approximation). The fluid 
is initially at rest, with the horizontal velocities Uz = Uy = 0 cm 
s“®. The physical conduction is calculated with an opacity fixed to 
Kc = 0.07 cm® g“® (see details in Cavecchi et al. 2013). 

On the horizontal components of the velocity we impose sym¬ 
metric boundary conditions at the vertical boundaries and reflective 
conditions at the horizontal boundaries®. For the horizontal compo¬ 
nents of the magnetic field we use antisymmetric boundary condi¬ 
tions both horizontally and vertically, while for the vertical com¬ 
ponent we use symmetric boundary conditions. By this choice, the 
field is forced to be perpendicular to the bottom and top surfaces. In 
all simulations presented here, we use horizontal and vertical reso¬ 
lutions of 192 and 90 grid points respectively. Our horizontal spatial 
resolution is therefore 3215 cm, as in Cavecchi et al. (2013), which 


^ We neglect electron screening. For comparison, a run which uses the full 
reaction rate including the screening factor (as in Fushiki & Lamb 1987) 
produces a flame velocity which is only 1.14 times faster. Since the focus 
of this paper is on the effects of the magnetic field we keep the simplified 
form Eq. (2). 

^ We ran some simulations without the damping force and obtained quali¬ 
tatively similar results, thus we conclude that this force does not affect the 
main conclusions of the paper. 

® Symmetric conditions imply that the values of the physical quantity are 
the same on either side of the boundary, while antisymmetric conditions im¬ 
ply that the values change sign. Reflective conditions imply that the veloc¬ 
ities have antisymmetric conditions while thermodynamical quantities like 
temperature and density are symmetric. In this way the boundaries on the 
left and right act like rigid walls that confine the fluid inside the simulation 
domain. 


is greater than the scale height. This illustrates the suitability of 
the hydrostatic scheme. Our numerical hyperdiffusive scheme has 
viscosity parameters vi = 0.03 and U 2 = 0.5 (see Braithwaite & 
Cavecchi 2012, for details). The corresponding values for the tem¬ 
perature and the helium fraction Y numerical diffusivities are taken 
to he 1 % of these®. We also use an explicit hyperdiffusive scheme 
for the magnetic field to ensure stability and avoid zig-zags, it is 
calculated as a multiple of the one used for the velocity (Braith¬ 
waite & Cavecchi 2012). 


2.2 Analytical calculations 


As shown in (Cavecchi et al. 2013), when the flame propagates in 
the absence of a magnetic field, the burning front is located along an 
inclined but nearly horizontal surface that separates as yet unbumt 
fuel from actively burning material (see e.g. Figure 1, left panel). 
A strong vertical shear is set up within the front. In an atmosphere 
without friction the horizontal velocity is in (quasi) geostrophic bal¬ 
ance with the horizontal pressure gradient. As a consequence, the 
vertical shear of the velocity is given by the thermal wind relation (a 
hurricane-like structure, see Pedlosky 1987 and SLU02 and Figure 
1 right panel): 

dvy ^ g dlogp 

dx f dz ^ ’ 

As the horizontal gradient of the density is concentrated within the 
burning front, so is the vertical shear of the horizontal fluid veloc¬ 
ity. When a magnetic field is present, this shear directly couples 
to the pre-existing vertical magnetic field Bx = Bq and generates 
the horizontal field By via ABy/dt = B^Ovy/dx; here d/dt is 
the Lagrangian derivative for the fluid element coming through the 
front. 

The horizontal field generates magnetic stress that back-reacts 
on the shear flow. The result is that the geostrophic balance is (par¬ 
tially) broken, which can lead to faster penetration of the hot burn¬ 
ing fluid in the top layers of the unburnt fuel. The acceleration of 
the front propagation due to the mechanical top-bottom momentum 
exchange was noted in SLU02 who introduced a phenomenological 
timescale ru for the frictional coupling and showed that the front 
velocity reaches maximum when rn ~ 1//. Here we study the ef¬ 
fect of the back-reaction of the flow-generated magnetic field and 
show that it does indeed accelerate the front propagation for a cer¬ 
tain range of magnetic field strengths (compare to sections 2.3 and 
3.7 of SLU02). 

Let Uh be the horizontal velocity at the top side of the burn¬ 
ing front. The vertical shear is, to order of magnitude, Vh/h where 
h is the vertical width of the front. The horizontal magnetic field 
generated in a fluid element upon its passage through the front is 

( 6 ) 

h 

Here Tn is the nuclear burning time which is also the characteristic 
time for a fluid element to cross the front. The horizontal magnetic 
field causes horizontal acceleration of the top velocity of the burnt- 
fuel 


BxBii 

~4n pH' 


(7) 


® Flowever, Runs 2, 3 and 7 needed a slight increase in the numerical ther¬ 
mal diffusivity to ensure numerical stability after some time during the sim¬ 
ulation, while Run 5 needed it since the beginning due to the fast fluid mo¬ 
tions (see Table 1 for the parameters of each run). 
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Run 

Bo (G) 

V (Hz) 

rR (s) 

rR (s) 

Rn (cm) 

(cm/s) 

Strength 

0 

— 

450 

1.77 X lO-'® 

— 

3.54 X 10® 

2.06 X 10® 

— 

1 

1 X 10’’ 

450 

1.77 X 10-® 

1.26 X 10“® 

7.91 X 10® 

7.28 X 10® 

Weak 

2 

1 X 10® 

450 

1.77 X 10--® 

1.26 X 10-® 

7.91 X 10® 

1.61 X 10® 

Intermediate 

3 

1 X 10® 

450 

1.77 X 10“® 

1.26 X 10“® 

7.91 X 10® 

4.78 X 10® 

Intermediate-Strong 

4 

1 X 10®° 

450 

1.77 X 10-® 

1.26 X 10-®" 

7.91 X 10® 

1.29 X 10® 

Strong 

5 

1 X 10'® 

11 

7.23 X 10“® 

1.26 X 10“® 

3.24 X 10® 

— 

Weak 

6 

1 X 10® 

11 

7.23 X 10-® 

1.26 X 10-® 

3.24 X 10® 

— 

Intermediate - Strong 

7 

1 X 10® 

11 

7.23 X 10“® 

1.26 X 10“® 

3.24 X 10® 

4.59 X 10® 

Strong 

8 

1 X 10®° 

11 

7.23 X 10“® 

1.26 X 10-®" 

3.24 X 10® 

1.33 X 10® 

Strong 


Table 1. Values of the initial purely vertical magnetic field Bq, the spin frequency u = Q/2 it, the rotational tr = 1/f and effective frictional (magnetic) 
timescales Eq. (14) defined in Section 2.2, the Rossby radius = y/gHr^ and the flame spread velocity ttf for the main runs discussed in the text. The 
last column reports the rough classification used in Section 2.3. The estimates for Tf^ are obtained substituting the value of Bq in Eq. (16), using = 0.1 

s and setting the numerical constants to 1. For we use g = 2 X 10^“^ cm s“^ and H = 10® cm. Vf is the value given by a linear fit to the horizontal 
position of maximum burning as a function of time. 


where p and H are the density and scale height of the burning fluid 
(Heng & Spitkovsky 2009). The overall horizontal acceleration of 
the burning fluid is caused by three physically distinct force fields: 
the horizontal pressure gradient, the Coriolis force, and the mag¬ 
netic back-reaction force: 


H ^ 

ah = aig — Cz 


fx Vh 


— <y.2 


Bq Tn ^ 

An pH h 


( 8 ) 


where ai and 0:2 are numerical coefficients of order 1, L is the 
horizontal extent of the front surface and we used — Bo, the 
seed field. To complete this equation, we note that first. 


ah OiQVh urn 7 


(9) 


and secondly. 


B ~ rourntt/i ■ €z' 


( 10 ) 


Here rbum ~ T^H/h is the timescale for the front to bum verti¬ 
cally through a patch of the ocean and ao is a constant of order 
1. Eq. (9) reflects the fact that the fluid has time Tbum/ao to ac¬ 
celerate from zero to Vh as the front bums vertically through the 
ocean, while Eq. (10) reflects the fact that it is the velocity compo¬ 
nent Vh ■ Cz that is responsible for the hot fluid covering the unbumt 
cold fluid. 

Eqs. (8), (9) and (10) form a closed system that can be solved 
for the two horizontal components of velocity Vh ■ Sz and Vh ■ Cy. 
The z-component is representative of the horizontal velocity of the 
burning front (compare SLU02, equation 34); it is given by 


V{ 


Vh - Bz = 


aigH 1/r 
. Tburn + I/t® 


( 11 ) 


Here 


r = 





ao 

"rburn 


-1 


( 12 ) 


where 


TA = H/va = ^/a^H/Bo (13) 


is the Alfven crossing timescale and va = Bo/2yTfp the Alfven 
speed. It can be seen that the effective frictional time is given in our 
case by 



a2'rburn 


(14) 


Let us consider several limits of this equation. In the non¬ 
magnetic limit, so that ta ^ Tbum, the front velocity becomes 
(using 1// < Tb urn) 


fnTT 

Vi ~ x/aiooy-, (15) 

J ^burn 

in agreement with Cavecchi et al. (2013) and SLU02. Note that vt 
is independent of the magnetic field coupling. However, for higher 
magnetic fields 


TA = 0.01 


H 


10® cm J \ 10® g/cm® 


1/2 


Bo 


10® G 


(16) 


may well be much smaller than rbum ~ few x 0.1s so that r ~ r^. 
So the next relevant comparison is that of the timescale rn with 
tr = 1/f. In the strong-field case where tr <C tr, which corre¬ 
sponds to 


Bo > Bo 


2.5 X 10® ( 5 ^ —- 

\ 10® cm / \ 10® g/cm® 

/ \l/2 Z a2rburn N~®/® „ 

V400 Hz/ V 0.1s / 


1/2 

X 

(17) 


the speed of the front becomes independent of the rotation rate: 


(18) 

V rbum rbum V OLO 

Einally, in the intermediate regime of Bo ~ Bo, i.e. rR ~ rR(^ 
rbum) the magnetic field strongly accelerates the speed of the flame 
propagation, reaching the maximum at Bq = Bo: 


max 


(wf) = y/ai/2 


x//rbum 


(19) 


The analytic arguments of this section are confirmed by the direct 
numerical simulations described below. 


2.3 Numerical simulations 

We simulated different configurations varying the initial vertical 
magnetic field intensity and the spin frequency of the star. We con¬ 
sider Bo in the range 10^ — 10^° G (broadly covering the range of 

^ Later, we will consider a case when v = 11 Hz. In this case 1// ~ 
7 X 10“® s, which still satisfies this condition. 
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magnetic fields in accretion powered X-ray pulsars Mukherjee et al. 
2015) and v = 450 Hz and v = 11 Hz; a summary of the parame¬ 
ters for each ran is reported in Table 1. We divide our simulations 
into three rough categories, which we call strong, intermediate and 
weak magnetic field based on the rn/rfr ratio and the develop¬ 
ment of the hurricane structure. We discuss these cases separately, 
comparing them to the case of Run 0, which has no magnetic field 
present and is rotating al v = 450 Hz. This simulation evolves 
as those discussed in Cavecchi et al. (2013). In particular, it has a 
well defined hurricane structure® identifiable by the thermal wind 
(see Figure 1, right panel). It also has a well defined burning region 
without turbulence (Figure 1, left panel). 

2.3.1 Strong magnetic field (10^° G) 

We start with the case of Bq = 10'^° G , Run 4. This field strength 
is at the upper end of the estimated burster magnetic field distribu¬ 
tion. In this case the initial overpressure at the top of the tempera¬ 
ture perturbation cannot bend the magnetic field lines significantly. 
Indeed, the frictional time scale is much faster than the rotational 
period (see Table 1). The initial motion is quickly suppressed and 
with no significant z velocity, no significant Coriolis force appears 
that can generate the thermal wind. 

In Figure 2 we show the results for Run 4. The top left panel 
shows the initial configuration of the magnetic field superimposed 
on the temperature profile. This configuration is the same for all 
simulations, apart from the value of the seed magnetic field Bq. On 
the right, we show the system at a later time: it can be seen that the 
field lines have not been bent significantly. The maximum bending 
is at the interface between hot and cold fluid, barely visible at z ^ 
3.6 X 10® cm. We observed the trace particle sliding down the hot- 
cold interface when the flame is approaching. When the flame has 
passed, they move mostly vertically along the straightened back 
field lines. Eventually, they only follow the expansion of the fluid, 
not moving anymore with respect to the fluid itself. 

We can compare this simulation to Run 0 (same spin, o = 450 
Hz, but without magnetic field) looking at the lower panels of Fig¬ 
ure 2 and at Figure 1. The first significant feature is that no hur¬ 
ricane structure has developed, as can be seen by the absence of 
any strong velocity Uy perpendicular to the plane in Figure 2. Ac¬ 
cordingly, the field has been bent in the y direction only slightly 
(By/Bx ~ 0.1 at its maximum). The only exception is a weak 
component of Uy right above the flame, where a component of at 
maximum 1.5 x 10® cm s“^ is seen briefly oscillating back and 
forth. We can see also a very weak component of the velocity at the 
top of the front, this is associated with the Coriolis force reaction to 
the motion of the fluid at the top of the interface. The second fea¬ 
ture is that after a steady state propagation has been established, the 
structure of the flame is very similar to the purely Coriolis confined 
one, even if the interface between hot and cold fluid is steeper, in 
the sense that there is no turbulence at the front, as opposed to what 
we see in the other simulations (see below). If we compare the es¬ 
timated values for ra and ra from Table 1 in the case of Run 4, we 
can see that the magnetic coupling should strongly dominate over 
the Coriolis force and therefore the length scale for the confinement 
should also be smaller than the Rossby radius. 

In Cavecchi et al. (2013) we discussed how the propagation 

® More precisely, the result of our simulations is an infinitely long roll, due 
to the use of flat Cartesian coordinates, which is the analogue to a hurricane 
structure. 


speed scales proportional to oc Rn/H, since i?R and H set the in¬ 
clination of the interface along which conduction can take place and 
the shallower the interface (longer Rb.) the more effective conduc¬ 
tion is. In the general case, Rb. should be replaced by the hot-cold 
fluid interface horizontal length scale, be it imposed by the Corio¬ 
lis force or the magnetic field or a combination of the two forces. 
As a confirmation, we find that the flame speed in the present case 
is slower than in the case without magnetic field since the flame is 
confined over a shorter distance. This simulation corresponds to the 
case of strong friction discussed in Section 2.2. The overpressure 
on the left is acting on all the layers of fluid ahead, due to strong 
magnetic vertical coupling, and the velocity is predominantly in the 
z direction, i.e. the horizontal direction of flame propagation in our 
notation. 


2.3.2 Weak magnetic field (llT G) 

We move now to the other extreme of Bq = 10^ G , Run 1. This 
approximately corresponds to the lower end of the magnetic field 
distribution for bursters. This case is much more dynamical. As can 
be seen in Figure 3 the late evolution of the flame and the field is 
strongly influenced by the motion of the fluid. Despite the fact that 
we call this case weak, it has to be borne in mind that even here the 
magnetic field is playing a role in the dynamics. 

At the beginning, the magnetic field is not capable of confining 
the initial pressure imbalance of the fluid and the fluid spills over 
sideways. At v = 450 Hz, the Coriolis force is strong enough to 
intervene and a hurricane structure begins to develop (Figure 3, 
upper left panel). The fluid motion drags the magnetic field along 
and strong horizontal components are generated. This agrees well 
with the fact that tr <C tr (see Table 1). However, the flame is 
very different from Run 0 (the case without magnetic field) as can 
be seen comparing the lower left panel of Figure 3 to Figure 1. 
In the case with weak magnetic field, we can see that the hot-cold 
fluid interface is much more extended. As a consequence, when we 
measure the propagation speed, we find that in the case of Bq = 
10^ G the flame moves faster than in the case of no magnetic field. 
Note that the value reported in Table 1 is the result of a linear fit to 
the horizontal position of maximum burning as a function of time. 
In this particular case, we found that the position of the flame was 
better fitted by a parabola than a line. This means that the velocity 
was accelerating till the flame reached the boundary. The value we 
report is therefore an average value. 

The reason for the longer interface and the higher speed is the 
reaction of the magnetic field which partly obstructs the Coriolis 
force. For example, it can be seen that the y component of the ve¬ 
locity, which characterizes the hurricane structure and ensures the 
Coriolis confinement, is smaller in absolute value than in the case 
with no magnetic field and is also limited to layers far above the 
level where most of the burning takes place®. The magnetic field 
component By is maximum at the flame location and steadily de¬ 
creases (apart from superimposed oscillations induced by waves, 
see below) once the flame has passed. The reduced confinement is 
responsible for the longer extent of the interface. This situation is 

® On the other hand, simulations with Bq = fO^ G , Bq = 5 X 10® G and 
Bq = 10® G , have velocities V{ = 2.06 X 10® cm s“®, Uf = 2.06 X 10® 
cm s“®, Vf = 2.14 X 10® cm s ^, which are almost identical to the one 
of Run 0. These simulations have a clear humcane stmcture without the 
presence of waves along the field lines. In these cases, the magnetic field is 
practically negligible. 
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Figure 1. Burning rate (left) and velocity in the y direction (perpendicular to the plane of the simulation - right) for the the non magnetic case, Run 0, with 
u = 450 Hz. The flame along the hot-cold fluid interface without turbulence and the thermal wind are cleai'ly visible. 


analogous to the case when friction is beginning to be significant, 
but it is still weaker than the Coriolis force. When the flame reaches 
the right boundary, the helium fraction is higher than in the case of 
no rotation or strong magnetic field. This can be easily understood 
if we consider that given the higher propagation speed, the hot fluid 
burnt for less time. 

Finally, we note that after the initial flame has developed, os¬ 
cillations set in along the hot-cold fluid interface, near the flame. 
Flere the baroclinicity is highest and Alfven waves are excited and 
move along the magnetic field lines. The waves propagate also 
backwards, at the base of the field lines. These, in turn, excite other 
waves along the field lines behind the flame, leading to the intricate 
configuration shown in the right, lower panel of Figure 3. 

There is an additional difference when comparing this case 
with the one of strong magnetic field or with the one of pure ro¬ 
tation; after ~ 0.3 s a peak forms in the height of the fluid (see 
Figure 3, upper right panel and lower panels). The position of the 
peak (z ~ 1.0 x 10® cm) does not seem to change appreciably 
for the rest of our simulation, and it is close to the position where 
our initial temperature perturbation was half of its maximum value 
(z ~ zo in Eq. 1). At this position the field lines are bent and re¬ 
main bent for the rest of the simulation. When the flame reaches 
the boundary on the right, the fluid is still burning and eventually 
the peak disappears. 

In order to confirm the correspondence with the initial pertur¬ 
bation, we ran a further simulation where the initial temperature 
distribution was given by Eq. (1), but with zo = 1.5 x 10® cm. 
We verified that also in this case a peak develops and it is now 
at z ~ 1.5 X 10® cm, i.e. where the centre of the new perturba¬ 
tion is. In both cases the magnetic field keeps memory of the ini¬ 
tial perturbation. The damping factor described in Section 2.1 is 
not the origin of the peak, since simulations with higher vertical 
resolution, the same physical parameters, but without the damping 
showed the same behaviour. The peak was there and even more 
pronounced. Inspection of the temperature distribution shows that 
the base of the fluid, below the peak, has a higher temperature (see 
Eigure 3, lower right panel). It appears that the heat is deposited 
in this region, mainly by viscous heating, following the waves that 


are localized below the peak. Our viscosity is mostly artificial and 
this may produce excessive heating, but if a similar heating process 
were proven to be physical, then such a peak could be the origin 
of an asymmetry in the emission during the burst rise and therefore 
lead to burst oscillations. 

2.3.3 Intermediate strength magnetic field (10® - 10® G) 

Einally, we address the case of intermediate strength magnetic field. 
The behaviour of Run 2, Bq = 10® G , is again different from the 
previous ones. In this case the initial heat reservoir perturbation that 
we impose on the temperature to trigger the burst is partly diluted 
during the very early stages of the simulation. The whole flame 
propagation is very chaotic. 

The partial dilution of the initial heat reservoir perturbation 
occurs in the following fashion. As can be seen in Figure 4, the hot 
fluid in the upper left comer of the simulation begins to slide to 
the right (panels A and B). The field lines follow the fluid, bending 
near the top at an approximate height of ~ 2 x 10® — 3 x 10® cm. 
Then, the curvature of the knee of the bent field lines increases and 
secondary, alternating, knees form below (panels B and C). The 
configuration becomes unstable around the upper knee (panel D) 
with the hot and cold fluid mixing: the field lines reconnect and the 
hot fluid moves right- and downwards, while the cold fluid goes up- 
and leftwards. This triggers a cascade of fluid rolls that propagate 
right- and downwards (panels D, E and F). Analogously, later in the 
simulation, we see a similar perturbation affecting the field lines. 
The perturbation propagates like a wave triggered where the weight 
of the fluid that is moving near the top of the simulation bends the 
field lines and these form a knee. The wave propagates from the 
formation of the knee down- and rightwards across the field lines: 
see the configuration of the field line around a ~ 1.8 x 10® cm 
in Figure 5 or the field line just after « ~ 2.4 x 10® cm at f ~ 
3.89 X 10~® s in Figure 6. This perturbation propagates towards the 
right of the simulation (see for example the field lines around 2 ~ 
4.8 X 10® cm in the second column of Figure 7), while additional 
waves are excited along the field lines that have been already hit by 
the front of the perturbation. Note that at t ~ l.llx 10~^ s another 
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Figure 2. Upper panels: initial and late conditions of the propagating flame in the case of Run 4, Bq = G and v = 450 Hz. The magnetic field lines 
are superimposed on the temperature profile. The magnetic field is strong enough to prevent most of the fluid motion. Even when the flame is propagating, the 
magnetic field lines are barely bent. Lower panels: nuclear reaction rate (left) and velocity in the y direction (right) for the same simulation. The absence of 
Uy shows that the propagation is driven by the sole magnetic confinement, since the Coriolis force is almost completely suppressed, to the point that a different 
color scale is needed to highlight the small components of the velocity. 


front of perturbation is now reaching 2 : ~ 4.8 x 10® cm (Figure 7), 
at f ~ 1.3 X 10“^ s it will be at 2 ~ 5.4 x 10® cm and it will reach 
the right boundary at f ~ 1.4 x 10“^ s. In order to find in which 
range of magnetic field strength the same initial bending of the field 
lines that leads to the beginning of the perturbation takes place, we 
ran a few more simulations and found that the same phenomenon is 
seen down to Bo ~ 3 x 10^ G, hut not anymore at Bo ~ 5 x 10® G. 
Note that this range encompasses the critical value Bo of Eq. (17). 

The partial dilution of the heat reservoir did not cause the 
flame to die, and the burning began right at the start of the simula¬ 
tion, but the flame is different from the well confined flame of the 
other cases we have treated so far. A pseudo steady state propaga¬ 
tion is established: from the main burning region, hot fluid tongues 
are launched up- and rightwards, these eventually turn downwards 
and leftwards when they reach regions where magnetic field lines 


have been compressed together. Secondary flames keep igniting 
where the front of the perturbation brings hotter fluid to regions of 
higher density. The new flame will propagate backwards till it joins 
the main burning region. First examples of these secondary burn¬ 
ing blobs can be seen in Figure 6. Clearer examples can be seen 
in Figure 7, since the burning is by now more intense also in the 
new flames, where the association of the secondary ignition sites 
with compressed field lines is more clear. The fluid keeps rolling 
between regions where the field lines have been compressed, mix¬ 
ing hot and cold fluid. An example can be spotted noting the de¬ 
formed shape of the field lines between 2 ~ 2.4 x 10® cm and 
2 ~ 3.0 X 10® cm at f ~ 1.17 x 10“^ s in the lowest right panel of 
Figure 7 or in the right, top panel of Figure 8 between 2 ~ 4.2 x 10® 
cm and 2 ~ 4.8 x 10® cm. 

These rolling eddies appear where the fluid is more baroclinic 
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Figure 3. Results for Run 1, Bq = 10^ G and u = 450 Hz. Upper panel: velocity perpendicular to the plane Uy at an early (left) and later (right) stage. 
Initially the height peak described in the text is not present: it develops at ~ 0.3 s, when the flame has passed the peak horizontal position. Lower panel: 
burning rate (left) and temperature profile (right) in the steady state propagation regime. Field lines are overplotted on the temperature profile. 


and driven by the baroclinic version of buoyancy. In this sense they 
are a form of convection (Pedlosky 1987). The fluid follows the 
compressed field lines on its left when going upwards, and the 
motion that follows after the fluid is stopped by the denser field 
lines ahead is reminiscent of the behaviour of the Parker instability 
(Parker 1966). This is understandable when considering the con¬ 
figuration of the field lines in this simulation. The field has been 
stretched and has developed a strong horizontal component. When 
the fluid at the top is stopped, it rests on these stretched lines, which 
are naturally subject to the instability. The fluid slides down and 
backwards along the inclination of the lines. When it stops, this 
fluid has higher entropy than the fluid above it and consequently 
is pushed vertically upwards again, generating the rolling motion 
between the regions of denser field lines. The continuous ignition 
of the fluid repeats this process until the right boundary is reached 
and the cycle cannot proceed anymore. 

The flame proceeds forward very fast, Vi ~ 1.6 x 10® cm s“^. 


in this fashion, until hy t ~ 3 x 10“^ sit reaches the boundary, 
where the field lines are being compressed. Finally, the fluid tem¬ 
perature reaches T ~ 10® K everywhere and by t ~ 5.98 x 10“^ s 
the fluid is expanding vertically approximately homogeneously and 
the field lines are straightening up. We note that these propagation 
timescales are in very good agreement with the observed rise times 
of observed bursts. 

Given the very fast initial phase where the flame reaches the 
extent of the boundary while the temperature does not rise above 
10® K, most of the fluid is still unbumt and it will be burned during 
the subsequent homogeneous expansion of the whole layer. Com¬ 
paring with the case of no magnetic field (Run 0) we find that dur¬ 
ing the whole propagation the hurricane structure is present, but 
the perpendicular velocity Uy is of smaller amplitude, (the maxi¬ 
mum absolute amplitude is roughly half the one of the case with 
Bo = 10^ G , Run 1, or of rotation only. Run 0). This is again in 
agreement with the fact that the flame velocity for Run 2 is faster 
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B, t = 3.67E-03S B, t = 8.49E-03s 








Figure 4. Run 2, u = 450 Hz and Bq = 10® G . The initial evolution of the thermal structure on the left added at i = 0 s to ignite the fluid. Shown are 
the magnetic field lines and tracer particles. The color of the particles is related to their initial horizontal position. The hot fluid that is trying to spill over 
is eventually blocked and directed downwards resulting in an early dilution of the original heat reservoir. The arrows indicate the horizontal position of the 
perturbation described in the text. 


than for Runs 0 and 1. The smaller the perpendicular velocity, the 
smaller the Coriolis force and the greater the propagation speed. In 
the left panel of Figure 9 we show the profile of the perpendicular 
velocity Uy, where the thermal wind is perturbed by the reaction of 
the Coriolis force to the horizontal 2 motion of the rolling eddies 
and the Alfven waves excited along the field lines. Finally we note 


that during the flame propagation after t ~ 8.46 x 10“^ s (Figure 
6) there is a temperature (and height) peak that is moving with the 
flame. The thermal wind is positive (out of the plane) at the flame 
front and negative (into the plane) behind. This propagation would 
appear as an expanding ring on the surface of the neutron star, at 
least until the effects of curvature become important. This is dif- 
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ferent from what we saw in the case of weak magnetic field, since 
there the peak was remaining at the same horizontal position, while 
the burning region expands, but also this configuration may lead to 
burst oscillations during the burst rise. 

We also ran a simulation with Bo = 10® G (Run 3, see Figure 
8 lower panels). Despite the fact that the simulation with Bo = 10® 
G is much more chaotic than this one and that no dilution of the ini¬ 
tial heat reservoir takes place when Bo = 10® G , these two cases 
share some similarities. The steady state configuration of Run 3 has 
a burning region which is more homogeneous than the one of Run 
2, but still it shows secondary burning bubbles and rolling eddies 
at the front (visible at 3.6 x 10®cm < z < 4.2 x 10®cm in the 
lower panels of Figure 8). This simulation bears similarities also 
with the case of strong field, since the spilling of the fluid is more 
controlled than in the cases of lower magnetic field. Flowever, the 
front is much more elongated than in that case, and flame propa¬ 
gation is consequently faster. The flame is faster also than in the 
case of no magnetic field. Run 0, but slower than in the case of 
Bo = 10^ G , Run 1, and Bo = 10® G , Run 2. The perpendicular 
velocity Uy is even smaller than in the case of Bo = 10® G (Run 
2), which is to be expected since the magnetic stress back-reaction 
increases with the seed magnetic field intensity. Moreover, the Uy 
velocity field at the front is oscillating and no real thermal wind is 
present (right panel of Figure 9). This is again very similar to what 
we see in the case of strong field, so that the case of Bo = 10® 
G should be regarded as a transition between the intermediate and 
strong field regimes. This behaviour is consistent with the fact that 
in this case m < tr (Table 1) which implies that at this point the 
magnetic coupling has become stronger than the Coriolis force. 

Run 2, Bo = 10® G, is the configuration that yields the fastest 
propagation speed among our configurations. This corresponds to a 
situation where the frictional coupling is most effective in increas¬ 
ing the speed of the flame and the speed value is near its peak. 
Indeed, Bo = 10® G is very close to our estimate of the critical 
field Bo. SLU02 discussed how in this case the flame propagates 
so fast that the ignition would be almost horizontally uniform: that 
is indeed what we observe. Increasing the field strength even more 
increases the coupling further, and the flame velocity will begin to 
decrease. The speed is expected to decrease because the pressure 
imbalance at the top has to push both the top and lower layers that 
are strongly coupled by the magnetic field. This is what we see hap¬ 
pening in the case of Run 3, Bq = 10® G and more drastically in 
the case of Run 4, Bq = 10^® G . 

2.3.4 Slow rotation - the case of IGR J17480 

Here we discuss the special case of a star spinning at t/ = 11 Hz, 
the rotation rate of the source IGR J17480 discussed in Cavecchi 
et al. (2011). Recall that in this case the Coriolis force is not capable 
of confining the fluid. We ran simulations using the same fields as 
we used for the case v = 450 Hz, so that comparison will be easier, 
but the concept of weak, intermediate and strong field should now 
be revisited: apart from the case with Bq = 10^ G , all fields have 
frictional time scales smaller than the rotational period and should 
be regarded as strong. 

First, we ran another simulation. Run 8, with a seed magnetic 
field of Bo = 10^® G , as in the case of Run 4. As expected, the ac¬ 
celeration due to magnetic stresses dominates that due to the Cori¬ 
olis force therefore controlling the flame propagation and the simu¬ 
lation results are almost indistinguishable from those at higher spin 
(Run 4). Therefore, also in this case, the flame propagation resem¬ 
bles that confined only by the Coriolis force at higher frequency 


(see Section 2.3.1) and thus we expect that a slow rotating source in 
this regime of magnetic field will still have a phenomenology sim¬ 
ilar to faster rotators without magnetic field^®. A simulation with 
Bq = 10® G , Run 7, also shows a behaviour very comparable to 
the case with the same field strength and o = 450 Hz, Run 3. That, 
again, can be understood in terms of the frictional time scale (see 
Table 1). Furthermore, we noticed that the thermal wind structure 
was absent u = 450 Hz, a fact that showed how the Coriolis 
force was no longer playing a significant role for the confinement 
and propagation. 

When the field is Bq = 10® G , Run 6, some differences 
between fast and slow rotation become more noticeable. After an 
initial phase where the fluid strongly bends the field lines, weak 
burning develops at the left end of the simulation, which triggers 
rolls that propagate towards the right similar to those of Run 2 
(Bo = 10® G , u = 450 Hz). These rolling eddies compress the 
fluid lines and climb over them at intervals. We see also fluid flow¬ 
ing down- and leftwards along the fluid lines when the eddies hit 
a new region of compressed lines. This propagation is similar to 
the case with Bo = 10® G and a = 450 Hz; however, at low ro¬ 
tation the length scale of these eddies is obviously longer (Figure 
10, upper panels). This velocity of propagation is very fast and at 
t ~ 2.6 X 10“^ s, the whole layer is horizontally approximately ho¬ 
mogeneously burning. It is difficult to track a neatly defined flame 
front, since secondary ignition sites ignite continuously due to the 
rolling eddies. Therefore, we do not report any velocity in Table 1, 
but a rough estimate gave vt ~ 3.5 x 10® cm s“^. Finally, when the 
burning is roughly homogeneous we see a convective pattern with 
three or four cells above the most intensely burning layer. These 
rolls disappear gradually during the rise of the burst, leaving at the 
end only two rolls, one per side, which eventually disappear as well, 
before the peak of the burst (Figure 10, lower panels). 

Finally, in the case of Bq = 10^ G , Run 5, we see that the 
fluid sloshes back and forth many times before ignition, which 
takes place at f ~ 1.94 s. This is to be expected, since nei¬ 
ther the Coriolis force nor the field (which is weaker) can signifi¬ 
cantly confine the fluid within the simulation domain and the waves 
bounce back only because of the reflecting boundaries. Without the 
magnetic field the fluid behaved very similarly, starting ignition 
at around t ~ 1.77 s. Convection develops above the horizontal 
layer that is burning, similarly to the case with Bq = 10® G and 
o — 11 Hz. Note that such late ignitions are probably caused by 
the finiteness of our simulation domain and the 2D dimensionality 
which prevents any dissipation in the y direction, while a neutron 
star spherical surface is much wider and in reality the flame would 
probably flume out. We conclude that the results of our simulations 
for the case of o = 11 Hz are in agreement with the predictions of 
Cavecchi et al. (2011): a field of Bq > 10® G is required in order 
to have flame confinement and propagation. 


3 SUMMARY AND CONCLUSIONS 

3.1 Summary of simulations at fast rotation (t/ = 450 Hz) 

In this paper we have presented 2D simulations of Type I bursts in 
the presence of both rotation and an initially vertical magnetic field. 

The velocity reported in Table 1 is slightly lower than the one for Run 
4, but that is just a consequence of the fact that we stopped the simulation 
early. When overplotting the position versus time of the two simulations, 
they agreed almost perfectly. 
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Figure 5. Continuation of the early evolution of Run 2 from Figure 4, Bq = 10® G = 450 Hz. Shown are the magnetic field lines with tracer particles 
superimposed. The arrows indicate the horizontal position of the later evolution of the perturbation described in the text. 


We have shown that in the case of a magnetic field of Bo = 10^° 
G , what we called the strong magnetic field case, flame propa¬ 
gation is independent of the rotation rates we used (u — 11 Hz, 
o — 450 Hz). The magnetic field prevents the development of 
any thermal wind (the hurricane structure described in SLU02 and 
Cavecchi et al. 2013) by inhibiting significant horizontal motions. 
The flame structure resembles very closely the one of a fast rotator 
in the sense that the fluid is not affected by turbulence, but propaga¬ 
tion is slower. On the other hand, we found that in the presence of a 
weak magnetic field (Bq = 10^ G ), the fluid confinement depends 
on the Coriolis force. In the case of fast rotation the presence of the 
field affects the structure of the front, making it more extended in 
the horizontal direction. This has the net effect of speeding up the 
flame propagation compared to a case with rotation only. 

A striking feature we found in the case of fast rotation with 
Bq = 10^ G is the development of a peak in the fluid height, due 
to high temperature, at a location approximately coincidental with 
the size of our initial heat reservoir in the temperature profile. A test 
with different size of the perturbation confirmed that the position of 
the peak is related to the initial conditions. From our simulations it 
looks like the origin of higher temperature is due to extra heating 
(mainly viscous) deposited by waves localized below the peak. Our 
viscosity is numerical, but if a similar mechanism were at work in 
nature, this would imply that there exist cases where a weak mag¬ 
netic field (when compared to the Coriolis force, in the sense de¬ 


fined in Section 2.3) would lead to the presence of an asymmetry 
in the emission pattern during the rise of the burst. This, in turn, 
could lead to the presence of burst oscillations. We also found that 
the velocity of the flame seemed to be accelerating until it reached 
the right boundary. Indeed, the flame propagation was better fitted 
by a parabola than a straight line. It would be interesting to check 
whether a constant velocity would ever be reached and on what 
length scale. In particular, it would be interesting to follow the evo¬ 
lution of this case on the spherical surface of a NS. 

The case of intermediate strength magnetic field with Bq = 
10® G at fast rotation showed a very remarkable series of results. 
First, we saw a temperature and height peak propagating with the 
flame, which would look like an expanding ring on the surface of 
the neutron star, until the effects of curvature would become sig¬ 
nificant. The duration and the possible effects on burst oscillations 
of this ring should be tested in 3D simulations. Secondly, we saw a 
perturbation triggered by the bending of the field lines at the very 
beginning of the simulation. However, this did not prevent igni¬ 
tion and flame propagation. Finally, this configuration showed that 
rolling eddies form at the flame front, which mix the hot and cold 
fluid and speed up the flame propagation by even an order of magni¬ 
tude when compared to the case without magnetic field. Also, sim¬ 
ilar eddies at the flame front form in the case of a field Bq = 10® 
G . In this case, however, there is no trace of the thermal wind and 
the hurricane structure that we still observe when Bo = 10® G and 
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Figure 6. Continuation of the early evolution of Run 2 from Figure 5, Bq = 10® G u = 450 Hz. Shown are the burning rate (left) and the magnetic field 
lines (right), with tracer particles superimposed. The arrow indicates the horizontal position of the perturbation described in the text at later time. 
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Figure 7. Continuation of the early evolution of Run 2 from Figure 6, Bq = 10® G u = 450 Hz. Shown are the burning rate (left) and the magnetic field 
lines (right), with tracer particles superimposed. The arrows indicate the horizontal position of the perturbation described in the text at later time. 
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Figure 8. The burning rate (left) and the temperature profile (right), with magnetic field lines superimposed of the later stages of the flame propagation for 
Run 2, Bq = 10® G V = 450 Hz, (upper panels) and Run 3, Bq = 10® G u = 450 Hz (lower panels). 


Bo = 10^ G . The burst proceeds with speed higher than in the 
case without magnetic field. 

The presence of the eddies at the flame front is very inter¬ 
esting, since they provide a means of heat conduction over a long 
length scale helping in speeding up flame propagation as compared 
to when only conduction is taking place (see for example the dis¬ 
cussion of Fryxell & Woosley 1982b, about the different speeds 
in the case of conduction, convection and turbulence). These giant 
rolls could also be responsible for bringing burning ashes to higher 
layers in the atmosphere, a fact that could lead to the generation of 
absorption lines (Weinberg et al. 2006). Their behaviour, in rela¬ 
tion with rotation rate and magnetic field strength, is very close to 
what is expected for convective rolls in the standard Benard prob¬ 
lem with both rotation and magnetic field (Chandrasekhar 1961). 
Indeed, their size decreases with increasing magnetic field (in our 
cases, going from Bq = 10® G to Bq = 10® G , for example) and 
increases with decreasing rotation frequency (as clearly seen for the 
case with Bq = 10® G going from u = 450 Hz to iz = 11 Hz). 


However, during the bursts we do not have a standard Benard con¬ 
figuration, since the fluid can be in motion due to the thermal wind 
structure (as in the case of Bq = 10® G ), the flame front is contin¬ 
uously moving and so is the source of heat that sets the temperature 
gradient. Furthermore, the configuration at the front is affected by 
the baroclinic motions, by the compression of the cold fluid under 
the weight of the expanding hot fluid and the Parker instability. The 
rolls clearly have an important role that could be amplified in 3D 
and further study is required. 


3.2 Summary of simulations at slow rotation (tz = 11 Hz) 
and implications for the source IGR J17480 

We performed simulations of a star with u = 11 Hz and magnetic 
fields in the same range as for the fast rotator case: 10^ 5; Bq 
10^® . Simulations with Bq = 10^® G and Bq = 10® G behaved 
practically identically to the cases with u = 450 Hz, since the mag¬ 
netic field was providing most of the confining force also in the case 
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Figure 9. The profile of the perpendicular velocity Uy in the case of Bp = 10® G (Run 2, left) and Bp = 10® G (Run 3, right) and u = 450 Hz. In the first 
case the thermal wind structure is still present, however disturbed by the reaction of the Coriolis force to the rolling eddy motions and by the Alfven waves 
along the field lines, while in the second case no hurricane has developed. Note that the color scale is different from the ones of Figures 1, 2 and 3. The white 
color patches correspond to velocity values above 5x10^ cm s“ ^. 


of fast rotation. The case with Bp = 10® G showed similar eddies 
to the analogue case at faster rotation, even though the length scale 
of the eddies was longer. The initial phase of flame propagation was 
even faster and eventually the explosion of the burst was mostly ho¬ 
mogeneously horizontal. In the case with Bp = 10^ G neither the 
magnetic tension nor the Coriolis force was strong enough to pro¬ 
vide any confinement and this led to the dispersion of the initial 
heat reservoir over the full horizontal extent of the simulations. We 
still saw a very delayed ignition, but we think that this is due to 
the smallness of our simulation domain compared to the size of the 
neutron star surface, a fact which prevents a more realistic dissipa¬ 
tion. 

We note that our simulations for the case of t/ = 11 Hz con¬ 
firm the conclusions of Cavecchi et al. (2011) about the case of 
IGR J17480. This source was discovered in 2010 in the globular 
cluster Terzan 5 (Bordas et al. 2010). It was found that IGR J17480 
displayed both accretion powered oscillations and burst oscillations 
at iz = 11 Hz (Strohmayer & Markwardt 2010; Altamirano et al. 
2010; Papitto et al. 2011). The accretion powered and the burst os¬ 
cillations had highly accurately matching phases and frequencies 
(Cavecchi et al. 201 1). This matching and the slow rotation rate led 
the authors to exclude the possibility that the Coriolis force could 
confine the flame and they also concluded that neither a Coriolis 
force confined hot-spot nor global modes could explain the burst 
oscillations. Based on approximated analytic calculations the au¬ 
thors estimated that the magnetic field of IGR J17480 could provide 
a dynamically effective confining force for the flame if B > 10® 
G. 

Therefore, a field weaker than Bp = 10® G would not confine 
the fluid and both ignition and burst oscillations should not be ex¬ 
pected. At Bp = 10® G , confinement is still weak and even if we 
do see the explosion leading to a burst, we do not expect to see any 
burst oscillation at this regime, since the flame propagates very fast 
and the whole domain ignites almost coincidentally. On the other 
hand, a field of Bp = 10® G or higher shows a confined propagat¬ 


ing flame, and we conclude that this is the minimum field intensity 
needed to have bursts that could display burst oscillations. 

3.3 Discussion of numerical and nuclear burning effects 

Before tackling the effects of the magnetic field, we will briefly 
discuss possible contributions that may affect the flame velocity 
we measure. 

One possible factor is resolution-related numerical effects. We 
could not perform full convergence tests for the simulations with 
magnetic field since that would have consumed many months, an 
amount of computational resources that we do not have at the mo¬ 
ment. However, the convergence tests in Cavecchi et al. (2013) 
showed that our code at similar resolutions to those used in this 
paper is converging, even if somewhat slowly. Furthermore, short 
convergence tests for the initial stages^^ of runs 0 — 4 showed the 
following trends. 

Run 0, being purely hydrodynamical, behaves exactly as the 
tests performed in Cavecchi et al. (2013) with the velocity of the 
flame decreasing slowly with resolution; run 4 behaves similarly, 
but the velocity at this stages seems already almost converged (ve¬ 
locity changed by only ~ 5 %). In the case of run 1 we note that 
at resolutions much lower than those in the text, many of the waves 
associated with the flame would not be resolved. The velocity mea¬ 
sured at this very initial stages is increasing with resolution, by 
~ 30 %, but the velocity we report in the text is evaluated at later 
stages that we did not reach in this short convergence tests. Runs 
2 and 3, on the other hand, should be more sensitive to resolution. 
This is linked to the fact that (a) the maximum flame position is dif¬ 
ficult to track for both simulations and (b) the convective rolls are 
captured differently at different resolutions. This situation is com¬ 
plicated by the fact that at different resolution our numerical diffu- 

^ ^ The extension of the tests much further that the initial stages is prevented 
by the limitations on computing time. 


MNRAS 000, 1-18 () 













16 Y. Cavecchi 


14.0 
12.0 
10.0 

E 8.0 

o 

C\J 

2 6.0 
4.0 
2.0 
0.0 

0.0 1.2 2.4 3.6 4.8 6.0 

10® cm 

B, t = 3.75E-01S 

14.0 
12.0 
10.0 

E 8.0 

o 

C\J 

2 6.0 
4.0 
2.0 
0.0 

0.0 1.2 2.4 3.6 4.8 6.0 

10® cm 




B, t = 2.00E-01s 




Figure 10. The evolution of Run 6, Bq = 10® G = 11 Hz, at different stages. Shown are magnetic field lines with tracer particles superimposed. Particles 
with the same color had the same horizontal position at the beginning of the simulation. Upper panels: the fluid is confined by the magnetic field only (left). A 
very fast initial propagation takes the form of rolling eddies (right). Lower panels: later stages of the burst when ignition is almost horizontally homogeneous. 
Convection sets in (left) and eventually subsides (right). 


sivities are different (Cavecchi et al. 2013). As expected, run 2 is 
more affected by resolution from the beginning, but the qualitative 
behaviour is the same as in our reference run, with the development 
of the initial perturbation described in Section 2.3.3. The velocity 
in run 3 seems almost converged in the early stages (difference in 
velocity is ~ 0.01 %), but we cannot exclude that it might change 
at later stages due to the onset of rolls (see Section 2.3.3). 

Other important factors are the physical burning and the opac¬ 
ity. Timmes (2000) performed a detailed analysis of ID horizon¬ 
tal laminar propagation of the flames, studying its dependence on 
the reaction network, thermodynamic conditions and conductivity. 
Timmes’ simulations were in good agreement with order of mag¬ 
nitude calculations predicting a velocity of the flame oc \/QJk, 
with Q the burning rate and k the opacity. A direct comparison to 
his results is impossible, mainly because the flame in our simula¬ 
tions, being vertically resolved, crosses different levels of pressure, 
while Timmes’ simulations were performed at constant pressure. 
Furthermore, we include cooling and fluid dynamics, which are ab¬ 
sent in that analysis. In Cavecchi et al. (2013) we proved that flame 
propagation in type I bursts can be expressed as a vertical flame 
velocity Vi± times a geometrical boosting factor set by the slope 
of the hot-cold fluid interface: ttf ~ V{±Rno/H. As we showed 
in the present paper, in many cases the net effect of the magnetic 


field is to interfere with the Coriolis force, therefore changing the 
slope of the interface and the boosting factor. Since this is the focus 
of the paper, we discuss it in the next section. However, a different 
conduction or a different burning regime as used by Timmes will 
change the vertical flame speed Our conduction is constant 
across the whole domain. While this is an approximation, it is less 
important than the reaction rate, since the value we use is calibrated 
to be the expected opacity at the conditions near the most burning 
regions. On the other hand, our reaction rate is limited to helium 
burning into carbon. A more extended reaction network with quick 
energetic reactions will boost the propagation speed by increasing 
the energy production rate at the vertical flame front. That is be¬ 
cause of the flame speed dependence oc v/Q (Timmes 2000). 
Slow reactions, like for example those limited by /(-decays in hy¬ 
drogen rich material, will probably have a smaller impact as long 
as there are faster ones occurring, since the slow ones will take 
place in the wake of the vertical flame front. Direct numerical sim¬ 
ulations are needed to better quantify these effects in the complex 
dynamical configurations of the type I bursts, in order to quantify 
the effects of composition and nuclear burning on the propagation 
speed and surface temperature distribution during all the stages of 
the bursts. This will also have a direct impact on the prediction and 
measurement of burst oscillations. 
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Figure 11. Left: velocity of the flame versus The blue diamonds represent the results for the fast rotator u = 450 Hz, while the red squares are the 

results for the slow rotator case (i^ = 11 Hz). The horizontal dashed line indicates the velocity of the flame when magnetic field is absent. Right: velocity of 
the flame versus Bq (note that Bg oc 1/t\ oc 1/The blue diamonds represent the results for the fast rotator t/ = 450 Hz and the red squares are the 
results for the slow rotator (i/ = 11 Hz) when tr Tfi. and the flame velocity is independent of rotation. The green circles are the results for the fast rotator 
when < Tfr and the velocity is affected by the Coriolis force. In these plots we include results both from Table 1 and footnote 9. 


3.4 The effect of the magnetic coupling on flame propagation 

speed 

The most remarkable result of our simulations is about the effect 
that the magnetic field has on the propagation speed of the flame. 
We summarize the results regarding the flame speed V{ in Figure 
11. Note that we include the values both from Table 1 and footnote 
9 (compare also to figure 2 of SLU02). 

In Section 2.2 we argued that there exists a field strength at 
which the propagation speed is maximal (Eq. 17, Bq ~ 2.8 x 10* 
G for u = 450 Hz). We also discussed how for decreasing field 
strength less effective coupling will lead to lower values of the 
flame speed, to the point that when tr <C tr the velocity should 
be independent of the magnetic field (Eq. 15). This is the trend we 
see in the left panel of Figure 11 looking at the results for u = 450 
Hz (blue diamonds). If Bg S> Bg, m 3> tr, the flame speed is 
expected to decrease continuously with the field strength. Again, 
this can be verified in Figure 11 . 

We also noted that when tr ^ tr the velocity of propaga¬ 
tion should depend only on the field strength via the frictional time 
scale TR (Eq. 18). This is more clearly seen in the right panel of 
Figure 11. There, the blue diamonds (v = 450 Hz) and the red 
squares (zz = 11 Hz) represents the results obtained in this regime. 
In particular, it can be seen that the velocities for Bg = 10® G and 
Bg = lO’^® G overlap almost perfectly. In contrast, for Bg ^ 10* 
G the results at zz = 450 Hz deviate from those at zz = 11 Hz 
because of the increasing importance of the Coriolis force on the 
formers (green circles), as expected. 

We conclude that our simulations confirm our analysis on the 
importance of the magnetic field in regulating the flame speed by 
providing an effective frictional coupling among the fluid layers 
and we note that for realistic field strengths we obtain speed values 
that are already in very good agreement with the velocities inferred 
from observations. 
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